Melting of two dimensional solids on disordered substrate 
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We study 2D solids with weak substrate disorder, using Coulomb gas renormalisation. The 
melting transition is found to be replaced by a sharp crossover between a high T liquid with thermally 
induced dislocations, and a low T glassy regime with disorder induced dislocations at scales larger 
than ^d which we compute (^^j ^ Re '^ Ra, the Larkin and translational correlation lengths). We 
discuss experimental consequences, reminiscent of melting, such as size effects in vortex flow and 
AC response in superconducting films. 



Although the phase diagram of the mixed state of high 
Tc superconductors iQ] has considerably evolved since 
their discovery, many questions remain. In d = 3 the 
low field region has been proposed as the topologically 
ordered Bragg glass phase [|| . The role of dislocations re- 
mains unclear near melting or at higher fields. To study 
analytically the effect of disorder on defective solids, 
d = 2 appears as a natural starting point, also important 
for numerous experimental systems besides thin super- 
conductor films |§-^, such as Wigner crystals in hetero- 
junctions m and on the surface of Helium M, magnetic 
bubbles arrays M. In the absence of disorder, the contin- 
uous melting of a 2d crystal occurs as dislocation pairs 
unbind at T°j and is well described by the KTNHY the- 
ory [ pOpl| ]. But no equivalent theoretical description of 
melting with substrate disorder exists, though it has been 
studied in many experiments 0,0,^E2l. 

Progress was made for the simpler problem of crystals 
with structural (i.e internal) disorder |13[ where melt- 
ing occurs at a temperature r,„ (a) shifted downwards by 
the disorder strength a (Fig.l). As discussed in [||, ||l^] 
is equivalent to treating only the long wavelength part 
a of the substrate disorder. On the other hand an ana- 
lytic RG study of the 2d Bragg glass including pinning 
disorder g, i.e short wavelength disorder (at the cost of 
excluding by hand dislocations) was performed recently 
M,15[. A complete treatment however should include 
dislocations. Indeed the general theoretical belief based 
on qualitative arguments and simulations [ p^|p7| ] is that 
no true solid (neither a lattice nor a vortex glass p^) ex- 
ist in 2D in presence of disorder at T > 0, and thus there 
should be no true melting transition. On the other hand 
signatures reminiscent of melting are observed in various 
experiments nM- This thus calls for further studies. 



In this Letter we derive renormalisation group (RG) 
equations describing a solid on a disordered substrate and 
allowing for dislocations. They generalize the KTNHY 
equations to weak pinning disorder near melting. They 
are obtained from the RG analysis, performed here for 
the first time, of the generic elastic vector electromag- 
netic coulomb gas (VECG) with disorder Pq |. Though 
the asymptotic RG flow is always towards strong coupling 
(disorder or/and dislocation fugacity) several studies are 
still possible for weak pinning disorder g but arbitrary a 



(see below) which strongly suggest that the pure melting 
transition is replaced by a sharp crossover between two 
very distinct regimes: (i) at higher T a weakly disordered 
liquid where disorder has no effect at large scale and dis- 
locations are present on all length scales greater than 
^° (T) as in the pure case (ii) at lower T a quasi Bragg 
glass regime, where disorder disrupts the quasi-long range 
order of the lattice beyond a length [|l^ Ra ^ Re and dis- 
locations appear on scales larger than a new length Cd{T). 
From the RG, here we estimate ^d(T) both for T < Tm 
where £,d ^ Re ^ a, and in the crossover region. At 
scales L < ^d, the low T defect free glassy regime is de- 
scribed by the theory of [|l^ , despite its defective nature 
in the thermodynamic limit. The sharp crossover pre- 
dicted between the two regimes at weak disorder is ana- 
lyzed and argued to account for experimental observation 
of 2D melting transition, and of the 2D peak effect. 
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FIG. 1. Schematic phase diagram: a and g are respectively 
the long and short (pinning) wavelength disorder. Inset: RG 
flow for weak g and T ^ T^ . The hatched region corresponds 
to the g = true solid which survives at <; = 0. 



We emphasize that our RG equations (5) are derived 
assuming that dislocations are thermalized. This as- 
sumption is reasonable near T,„ and allows us to ana- 
lyze the crossover near melting. A fully consistent RG 
analysis including pinning of dislocations which become 
important at lower temperatures O] , is beyond the scope 
of this paper EG]. However, as a first step, we also stud- 
ied the case of no pinning disorder 5 = at T < T^fj2 
by extending the analysis of (2^p2l on a simpler model. 
A further study of melting is presented elsewhere |lq| . 



In the absence of disorder and of dislocations a two di- 
mensional crystal is described by a smooth 2D displace- 
ment field u(r) and an elastic distortion energy: 



-He 



2ceeuf, + (cn - 2cQe)uiiUj 



(1) 



where cn and cee are respectively the compression and 
the shear elastic moduli, and Uij = -^{diUj -t- djUi). For 
simplicity we study a triangular lattice of spacing a. 
Next, disorder is modeled by a gaussian random potential 
V{y) with correlator V{y)V{t') = h{r — r') of range r/, 
and couples to the density of vortices p(r) — J^i '^('' ~ ^i) 
according to Tip = J^p{r)V{r). Decomposition of p(r) 
and of V(r) in Fourier components gives several couplings 
between u(r) and the disorder M. The two main contri- 
butions [ p3[ are a random stress field aij arising from 
the long wavelength part of the disorder, and a 'random 
phase field' c/ji, which comes from the part of the (pin- 
ning) disorder with almost the periodicity of the lattice: 



each nia belonging to the reciprocal lattice, with the con- 
straint jlj] ^^ nia = as well as the usual neutral- 
ity condition. Integration over the replicated field ug(r) 
leads to an elastic VECG (instead of the purely elec- 
tric studied in p4|) which in its most general form is 
defined hjZ^= E{,,„.} C'^ DLi / ^Y[h^,m^]e-^ 
where the integrals are restricted by \ra — i"/?! > a, the 
Y[h, m] are the composite charges fugacities and C{b, m) 
is a combinatorial factor. The action is given by 



S[b, m] = i<{ b'' * V^ini, n^) * b'^ 



(3) 



Tip 



= I (JiiU 



V "y 



2W^ J2 cos(K,.u(r) + 0,(r)) (2) 

i/=1.2,3 



with (e^('^-('-)-'^^'('-'))) = Si,^^,S^{r ~ r'), the K^ are the 
first reciprocal lattice vectors (of length Kq — 47r/V3a) 
and the correlator (Jij{r)aki{r') is parametrised Q by 
All, Age whose bare values are An — pghq^Q, Age — 0, 
while g — a? p^hq=Ka/T^ where po is the mean density. 
Note that if rf > a, g can be greatly reduced (e.g. by 
a factor exp(— c(r'//a)^)) with respect to An^ge- For 
out of plane disorder (such as in B), varying the rela- 
tive strength of the two types of disorder can be done by 
changing the distance between the lattice and the disor- 
der plane. Besides this contribution, an underlying dis- 
order potential induces local prefered orientation of the 
lattice pq l : this new random field (random torque) cou- 
ples to the local bond angle [^ 6 — \{dxUy — dyUx) ■ 

Tis = /r A{Y).e{T) with A{y)A{v') = 4A^(5r^r'- Evcn if A^ 
is null in the bare model, it will be renormalised to finite 
value (see (|5|)), and must be taken into account from the 
beginning since it is a new independent disorder strength, 
whereas without dislocations [MJ it can be absorbed in a 
redefinition of Agg (Agg -^ Agg -f A^). 

To describe plastic distortions of the lattice one spfits 
u = Uq -|- Ud into a smooth phonon part Uq and a dis- 
placement field Ud,i{r) = ■^Jj.,Gij{r — r')bj{r') due to 
edge dislocations of Burger vectors density b(r), where 

g,,{r) = S,,^r) + ^6,,G (-) + ^"~^^^ ,fcg..fc(r) 
Cn \a/ cii 

with G{r) + i$(r) — \n(x + iy) where r = (x, y). The in- 
teraction Hik{v) is defined by [|6| H.,k{r) = ^ - \5ik. 
Introducing n replicas and averaging over disorder using 
the Villain form for the cosine coupling in (H) we obtain 
| p^ , p^ in each site repficated charges m = (mi,..m„), 



denoting b*y*b = Ey/r^^,b,(r)y„(r-r')bj(r'). The 
interactions are ttV°-'^{ki,K2) = n\'^5ijG — n'^'^Hij and 
2TTW^f{n5,KQ) = %,5"^$ + ^re^jG + <%fei?^fc. The 
renormalization of this elastic VECG goes beyond the 
previous analysis of the EGG [g^ j due to the presence of 
new marginal operators corresponding to the elastic in- 
teractions y, W . Note that usual electric/magnetic self- 
duality for (g) now reads: ki <-^ K3 ; K2 <~^ ^^4 ; K5 ^^ 
— K5 ; Kg <-> Kg and exchanges direct and reciprocal lat- 
tices. The full RG equations and applications to various 
2D models (depending on the definitions of the Ki.g) 
is presented elsewhere Q. Here we study the model 
(0,0) for which the Ki...g are replica matrices of the form 
Kf^ = 7t7(5°'^ — Ai with the matrix definitions : 
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cqqC^i 


(4c) 



where the original parameters of the model (|l|,g) have 



'11,66 — 

yA^. Without 



been embedded in three replica matrices C 

011,66-^°" - yAii^gg and 7"^ = -fS'"' 

coupfing to a periodic lattice [g^, one sets 7 = 0. 

This VECG can be studied by the RG as in ^ by 
incrementing the hard-core cutoff a — *■ ae' in a three 
steps coarse-graining: reparametrisation of the inter- 
action, which gives the dimension of the operators : 
diy = (2 - g(Ki - Ai))y,dig = (2 - -g 2^3)3; fusion 
of electromagnetic charges separated from a < d < ae'; 
annihilation of dipoles of diameter a < d < ae\ which 
defines scale dependent interactions Ki...g(^) |1^. Here 
we restrict ourselves to charges of minimal fugacity: (i) 
charges b with single non zero replica component ba = G 
(a lattice vector of minimal length), y(b,0) = y (ii) 
charges m with two (opposite) non zero replica compo- 
nent nia — — rtic = Kj,, Y{0, m) — g pO[ . We show after 
rather tedious calculations |jlq| that (Kl) is renormalisable 
to lowest order within the form (0), needed for consis- 
tency. We obtain the general RG equations, which, in 
the case of (QJ^), generalize the KTHNY equations: 



dig = 



'2-—(—- ^ + '^^^ 



4tt C66 



4c, 



'66, 



g- Bm{a)g^ 






(5a) 
(5b) 

(5c) 
(5d) 



g^((4c66-if)^(/o + /i)+i^%).g^ (5e) 

where /o,i = /o,i(a) and /o,i = /o,i(a) are modified 
Bessel function. We liave defined K — 4Tki — 4Tk2 — 
4(cii - C66)c66/cii, aKy2 = T2(Ai + A2)_= A66(l 
2C66/Cii) + Aii(c66/cii)2, ,5 = 4A, 



^, and yjl B„i(a) = 
2^(2/o(a/2) - /o(a)), a = f^(if/r - aK^T^ + S/T^) 
and a = p^K^KT /{1&ttc\q). We introduce an additional 
parameter p analogous to the p-fold symmetry breaking 
field in |||, setting Kf^ -> pK^^ in (|). The physical 
model corresponds to p = 1. We emphasize again that 
(h) should be valid at high enough T (near melting). 

We start by the simpler case of no pinning disorder 
g = 0, where further extensions of (ra) at low T can also 
be given. The modifications induced by g > are dis- 
cussed later for T > Tm/2. g ~ is experimentally 
relevant when disorder (An^gg^-y) varies very smoothly 
at the scale of the lattice. At high enough temperature 
T > Tm /2, eqs. (|^) show that the solid is stable at weak 

(T + 5/K'^ and at low tem- 



75^^ 2 



disorder a < tA— where aa 

647r 



peratureT < T^{a,d) = -^§1^(1-^1-6477^) (note that 

All, 66, 7 are unrenormalized) . At T — Tm{(T,S) < Tm it 
undergoes a true KT like melting transition (see Fig.l) 
to a high temperature phase where dislocations prolif- 
erate. The correlation length at the transition is given 
by C+ ^ exp(const(T — Tm)~'') where v, computed in 
|18| depends continuously on a and S and vanishes for 
W = g^. When the random shear and torque are null 
Aj = Age = 0, one recovers Nelson's results [Q. At 
lower temperatures (T < T^/2), freezing of dislocations 
leads to modified RG equations for T = T/{2o?aK) < 1. 
We showed |1^ that it is possible to extend to the elastic 
model the approach of 123] for the simpler scalar model, 
by defining the appropriate dislocation fugacities y : 



diy^{2- 



1 



y , diS^O 



32ttW^ 
di{TK-^)^C{l-q)f , dia^Cqf 



(6a) 
(6b) 



where g w T is the ratio of frozen dislocations [^ . The 
solid phase thus survives for an < q^ , as in [22] , and the 
physics is dominated by rare favorable regions for frozen 
dislocations |Q, leading to the phase diagram of Fig. 1. 
Positional correlations decay algebraically with exponent 

_ GfT/ -1,-1 ,Ai 
— 4ii- >^'-ll +''66 +T^ 



VG 



^ ,y ^ ) , which leads at the 
transition T ^ T^/2 to ^ < rjc < ^ (for Age = A.^ = 0), 



We now turn to the effect of pinning disorder. Let us 
first recall that in situations where dislocations can be 
neglected (see below) one sets y ~ in (|^) and recovers 
the Bragg glass phase [Q g*B™(a) = 2(1 - T/Tg) = 2r 



which exists for T < T„ 



., - «^(2 - KViK,))-' 
(much larger than Tm)- Due to the unbounded increase 
of All gg ~ ^: disorder induced displacements grow as 
u ~ Inr. The rest of the paper is devoted to study- 
ing the situation where both weak pinning disorder and 
dislocations are included, using eqs.(||). The RG flow 
shows a sharp crossover between two distincts regimes 
(see fig. 2) : a high temperature regime characterized by 
a correlation length S,\{T) unaffected by pinning disor- 
der {i.e the same as given above), and a low temperature 
glassy regime where translational order decays beyond 
Ra ^ Re, though equilibrium dislocations are separated 
(apart from small dipoles of size ~ a) by the larger length 
Cd(T) (defined by Cd = ae^\y{l*) - 1). We also study 
the crossover by parametrizing temperature and disorder 
using the correlation lengths f ^ (T) of melting in the ab- 
sence of pinning disorder and the Larkin length in the 
absence of dislocations R^ ^ g^~ (see (35) in Q). 
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FIG. 2. The length ^d/a, as a function of temperature, for 
various pinning disorder strength g, and (rescaled) scale de- 
pendant parameters in the two regimes (inset). 



In the low T glassy regime as I increases y{l) starts by 
dropping quickly to a negligible value (Fig. 2) whereas g{l) 
increases and reaches for L = ae' ^ Re ^ R^ a. plateau 
value corresponding to the Bragg glass fixed point g* de- 
fined above [|14|J30[ . Within this range of scale K and 
C66 remain constant since y(l) ~ 0. However the steady 
growth of a and 5 {g{l) being finite) eventually makes the 
eigenvalue of y positive. Since y{l) is then very small, 
it takes a large I* before it reaches ~ 1 and explodes 
to a non perturbative value. Thus dislocations, which 
are strongly suppressed in the intermediate Bragg glass 
regime, appear only on very large length scales: to esti- 
mate S,d, we can neglect the y^ terms in (H) {y{l) ^ 1 up 
to ^rf). Integrating the corresponding RG eqs., we obtain 
in the limit of weak disorder (large i?°) 



u 



R^e"^^^' 



with b = 1t\/ \{(Jc — 



c^oj 



(7) 



pil where A is a constant and Re ^ it!°(c66(^c)/c66(0))~ 
is the true Larkin length. Note that Re ~ i?" i'^ ^^i^ ^^'^ 
T regime, except very near the crossover region. 

In the high T liquid regime, y{l) increases quickly, caus- 
ing g{l) to decrease to 0. Thus the correlation length is 
the same as without pinning disorder S,\. However this 
regime cannot be considerer as an hexatic phase: since 
the renormalised A-y is nonzero (computable from our 
RG), the corresponding random field destroys |2^ long- 
range orientational order |18| , p2[ beyond Rh ^ [R^^Y /a 
(since A^ ^ g^) leading to a liquid at large scales. 

Finally, the crossover between the two regimes is dom- 
inated at weak enough disorder by the 5 = fixed point 
(see fig.l). Studying the RG flow around the KTNHY 
separatrix [|ll[ we find ^d va the crossover region mm : 



id 



R° 



{InRoy^ 



It 



Inflg 
In^^ 



(8) 



where the ± sign depends on T > T^ or T < T,o [|| and 
T w 0.8, v as in |1J] near pure melting (ct « 0). When 
increasing pinning disorder (or decreasing temperature), 
^d gradually goes from (||) to the low T behaviour (^. 

This sharp crossover will have consequences for finite 
size effects in 2D systems such as thin superconducting 
films P,plH]. In experiments of narrow channel vortex 
fiow in NbsGca films H, one probes the visco-elastic re- 
sponse of the lattice for T < T„i on scales much smaller 
than ^d (see fig. 2) : the system responds as a solid with 
a finite shear modulus cge, and one observes a "diverg- 
ing" correlation length (or viscosity H) when approach- 
ing TmiT), as without pinning. In much larger system 
for T < Tm, the response of the lattice should in fact be 
liquid- like with a large viscosity v ~ Cd(^)- I^ ^^ ^^" 
perimen ts ]g] , the 2D vortex lattice is probed on a length 
leu ~ \/D/uj. By varying to, a crossover at l^i^ = R^^ 
was observed at low T in l^l between a low oj liquid- 
like behaviour and an activated glassy behaviour. The 
length R^^ strongly differs [|6| from the corresponding Re 
extracted from critical current experiments. We argue 
that the observed R,^^ of Q corresponds to the correla- 
tion length £,d{T) defined here. Indeed the scaling of R, 
with the sample thickness d (i.e disorder strength) in 
is consistent with (^. Finally, using (||) we compute the 
scale dependent ceeil) and obtain the softening of cee on 
scale Re by dislocations. The self consistent definition of 
Re given above then leads to a quantitative description 
of the increase of critical current (peak effect) in 2D films 
| fi2[ near melting, further studied in p8[ . 

To conclude we extended the KTNHY analysis in pres- 
ence of weak pinning disorder, predicted and analyzed a 
sharp crossover near pure melting. To go beyond our 
study in the low T region would necessitate a controlled 
RG method pi] to describe frozen topological defects. 
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